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Abstract 

We consider a stochastic model, describing the growth of two com¬ 
peting infections on The growth takes place by way of spherical 
outbursts in the infected region, an outburst in the type 1 (2) infected 
region causing all previously uninfected points within a stochastic dis¬ 
tance from the outburst location to be type 1 (2) infected. The main 
result is that, if the infection types have the same intensity, then there 
is a strictly positive probability that both infection types grow un¬ 
boundedly. 
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1 Introduction 

In Deijfen (2003), a growth model is introduced that describes the spread of 
an entity, henceforth referred to as an infection, on W^. The model is gener¬ 
alized in Deijfen et al (2003) to encompass two infection types, making it a 
model for competition on Deijfen et al (2003) is mainly concerned with 
the case when the infections have unequal intensities and the main result is 
roughly that two infections with different intensities cannot simultaneously 
grow to occupy infinite parts of In this paper we consider infections with 
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the same intensity and show that in this case there is a positive probability 
that both infection types reach points arbitrarily far from the origin. 

The Richardson model, introduced in Richardson (1973), describes growth 
on Sites are either healthy or infected, the dynamics being that a healthy 
site is infected at a rate proportional to the number of infected neighbors 
and once infected it remains so forever. The model of Deijfen (2003) can 
be viewed as a kind of continuum analog of the Richardson model. Starting 
from an arbitrary initial set, the growth is driven by randomly sized spherical 
outbursts in the infected region, whose value at time t is denoted by St- 
The time until an outburst takes place is exponentially distributed with 
parameter proportional to \St\, where | • | denotes Lebesgue measure, and 
the outburst point is uniformly distributed over St- When an outburst 
occurs it causes the previously uninfected parts of a ball centered at the 
outburst point to be infected, the radii in the outburst balls being i.i.d. 
random variables with distribution F. The main result for the model is a 
shape theorem stating that, if 

e~‘^^dF{r) < oo for some (f <0, (1) 

then, on the scale 1/t, the set St has an asymptotic shape, which by rota¬ 
tional invariance is forced to be a Euclidean ball; see Theorem 1.1 in Deijfen 
et al (2003). This is a stronger result than for the Richardson model, for 
which the existence of an asymptotic shape is established but the nature of 
the shape so far remains unknown. 

In Haggstrom and Pemantle (1998,2000) a generalized version of the Richard¬ 
son model is introduced that contains two infection types, referred to as type 
1 and type 2 infection respectively. The dynamics is that, for i = 1,2, a 
healthy site becomes type i infected at rate A* times the number of type i 
infected nearest neighbors and once infected it never recovers. Here A* G M 
indicates the intensity of the type i infection. Furthermore, when a site is 
struck by one of the infection types, it is assumed to become immune to the 
other type, making the process a model for competing growth on Z'^. 

The continuum model of Deijfen (2003) can also be extended to incorporate 
two infection types; this is done in Deijfen et al (2003). Like in the one-type 
model, the growth takes place by way of spherical outbursts, the type of 
an outburst being determined by the region in which it takes place. Let 
St (St) denote the region occupied by type 1 (2) infection at time t. The 
time until an outburst occurs in SI {i = 1,2) is exponentially distributed 
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with parameter Ails'll, where A, € M is the intensity of the type i infection, 
and the outburst point is chosen uniformly in SI- An outburst in the type i 
infected region hurls out type i infection to all previously uninfected points 
in a ball centered at the outburst point. The radii in the outburst balls are 
i.i.d. random variables with the same distribution, denoted by F, for both 
outburst types. Note that the infection types are assumed to be mutually 
exclusive, that is, a point can be contained in at most one of the sets and 

Si 

According to the above mechanisms, the two-type models describe a com¬ 
petition for space - on and respectively - between two infections. For 
both models there are three possible results of this competition; Either one 
of the infection types wins or there is a draw. For i = 1, 2, define 

Ai = {the type i infection reaches sites arbitrarily far from the origin}. 

We say that the type 1 (2) infection has won if - 1 A 2 (“'Ai) occurs and with 
a draw we mean that A = Ai n A 2 occurs, that is, a power balance is 
established so that both infection types keep growing indefinitely; here -■ 
denotes set complement. 

In the discrete model and in the continuum model with bounded outburst 
radii, clearly both events -lAi and - 1^2 - and thereby also the events Ai 
and A 2 - have positive probability regardless of the values of the intensi¬ 
ties Ai and A 2 . In the continuum model with unbounded support for the 
outburst radii, this fact requires a proof; see Proposition 5.1 in Deijfen et 
al (2003). The mathematical challenge with the two-type models lies in 
deciding whether the event A has positive probability or not. Intuitively, 
the state of equilibrium represented by A should be possible to maintain in 
the long run if and only if the infection types have the same intensity. To 
state this conjecture more formally, let denote the probability law for 

a discrete two-type process started from two single nodes, one at the origin 
and one at the point (1, 0,... , 0), and write A 2 fQj- distribution of a 
continuum two-type process started from two unit balls located next to each 
other, one centered at the origin and the other at the point (2,0,..., 0). 

Conjecture 1.1 For any dimension d>2, we have 

(a) >o^Ai = A 2 ; 

(b) if F satisfies then P^^’^^{A) > 0 Ai = A 2 . 
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Remark 1.1 It is shown in Deijfen and Haggstrom (2003) that the possi¬ 
bility of infinite coexistence in the lattice model is not affected by the initial 
configuration, that is, whether is positive or not does not depend 

on the initial state of the model, given of course that neither of the infec¬ 
tion types surrounds the other already at the beginning. In Deijfen et al 
(2003) an analogous result (stated as Proposition 13.11 in the present paper) 
is proved for the continuum model provided the technical condition 

F(e) > 0 for all e > 0 . ( 2 ) 

Hence the particular choice of initial sets is unimportant in Conjecture 1.1. 
□ 

Efforts have been made to prove Conjecture 1.1 and a number of partial re¬ 
sults have been obtained. As for the case with unequal intensities, it is shown 
in Haggstrom and Pemantle (2000) that, if Ai is fixed, then = 0 for 

all but at most countably many values of A 2 , and in Deijfen et al (2003) the 
same result is established for pPA 2 under the assumption that ([T]) holds for 
F. Unfortunately, to strengthen these results to cover all intensities A 2 7 ^ Ai 
seems to be marred by large technical difficulties. 

The discrete model with equal intensities was first studied by Haggstrom 
and Pemantle (1998). Their main result is that, if Ai = A 2 and d = 2, 
then P^^’^^{A) > 0 , that is, on a two-dimensional square lattice, mutual un¬ 
bounded growth is possible for infections with the same intensity. Recently 
Caret and Marchand (2004) were able to extend this result to arbitrary di¬ 
mensions. While the technique of Haggstrom and Pemantle relies on some 
very specific features of the square lattice the proof of Caret and Marc¬ 
hand is much less sensitive to the exact lattice structure, and thus more 
amenable to generalizations. In this paper, we exploit the ideas of Caret 
and Marchand to establish infinite coexistence in the continuum setting. 
That is, we prove the following. 

Theorem 1.1 Let d > 2 and assume that F satisfies conditions m and 
Then, if Xi = X 2 , we have P^^'^'^{A) > 0. 

Remark 1.2 Condition ([2]), which says that arbitrarily small outbursts hap¬ 
pen with positive probability, is only needed in order to apply Proposition 
l3.1l below (quoted from Deijfen et al (2003)). Without assuming ([2]), we still 
have the slightly weaker result that infinite coexistence has positive prob¬ 
ability for certain other compactly supported initial states - see Theorem 
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[Q Furthermore, we strongly believe that Proposition 13.11 holds without 
the condition ([2]). □ 


The most important step we need to take in order to be able to invoke the 
Garet-Marchand ideas is a new coupling construction of the two-type con¬ 
tinuum model, designed in such a way that certain subadditivity properties 
of the process are reflected in the coupling. This construction is described in 
Section 2 together with some preliminary results, and the proof of Theorem 
II. Il ls given in Section 3. Superficially, our proof may not seem to have much 
in common with the Garet-Marchand paper, but the underlying key idea is 
the same. Translated to the continuum setting, the basic idea is to observe 
(|7|l and to demonstrate that (A) = 0 would lead to consequences that 

contradict dZD- 

2 A new coupling construction 

Let us begin by describing in more detail how the continuum two-type model 
with equally powerful infections can be generated. To be able to exploit the 
Garet-Marchand ideas, we need to modify the construction in Deijfen (2003) 
and Deiifen et al (2003) in order to make a certain subbadditivity relation 
(HD hold. 

To construct the model in d dimensions, consider a Poisson process on 
with rate A, where the extra dimension represents time and A denotes the 
common intensity of the infection types. Associate independently to each 
Poisson point a random radius variable with distribution F and pick two 
bounded disjoint regions Sq,Sq C to start up the growth, the set Sq 
( i = 1,2) indicating the initially type i infected region. 

Starting at time 0, the growth is now brought about by scanning within the 
region (S’)) US'q) x M upwards along the time axis until a point in the Poisson 
process is hit. An outburst then occurs at this point, causing the previously 
uninfected parts of a ball around the outburst point to be infected. The 
radius of the outburst ball is given by the radius variable associated with 
the Poisson point and the type of the transmitted infection is determined by 
in which one of the regions Sq x M and S'q x M the Poisson point is found. Let 
Ti denote the time at which this first outburst takes place and, for z = 1, 2, 
let denote the type i infected region after the outburst. Here the region 
with the same type as the outburst might be enlarged as compared to before 
the outburst, while the other region remains unchanged. 
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To generate the next outburst, we keep scanning upwards in time for Poisson 
points, but in such a way that the region Sq U Sq and the region U ) \ 
(S’)) U S'q) are treated somewhat differently: for the region Sq U S'q, we scan 
from time Ti and upwards, while simultaneously for U \ (<50 
we scan from time 0 and upwards. Hence, the time T 2 of the next outburst 
is given by the minimum of T 2 and T 2 , where T 2 is the time coordinate of 
the second Poisson point in (S'q U Sq) x M, and T2 is the sum of Ti and the 
time coordinate of the hrst Poisson point in U \ (Sq U Sq) x M. 

Subsequent outbursts are generated by similarly scanning up the time axis. 
This means that a Poisson point at coordinates {x, t) with x G and t > 0 
represents an outburst at x, not at (absolute) time t, but precisely t time 
units after the point x first became infected. We thus obtain two increasing 
sequences {<5(\j)} and {<5(^„)} of regions in indicating the type 1 and type 
2 infected regions after n outbursts, and a strictly increasing sequence {T„} 
specifying the time points for the outbursts. For t € [rn,rn+i), the type i 
infected region at time t is given by SI = and the total infected region 
at time t is St = S} U Sf. Clearly this dehnes a model where the time until 
an outburst occurs in Sj is exponentially distributed with parameter A|5j| 
and the location of the outburst is uniformly distributed in SI, as desired. 

In summary, the probability space underlying our construction, is simply a 
marked Poisson process on with rate A and i.i.d. marks with distribu¬ 
tion F. A more thorough description of the two-type model can be found 
in Section 3 in Deijfen et al (2003). The construction there is different from 
the one given here, but it is easy to see that, if the infection types have 
the same intensity, the two construction methods result in growth processes 
with the same distribution. The difference between the two constructions 
matter only when we couple realizations of processes starting from different 
initial conhgurations. 

If we do not distinguish between the type 1 and the type 2 infection, then 
the symmetric two-type model reduces to the one-type model of Deijfen 
(2003), that is, the total infected region in the two-type model behaves like 
a one-type process. An important consequence of this is the existence of a so 
called time constant, denoted by /r, indicating the inverse asymptotic speed 
of the growth. Let B{x,r) denote the closed ball with radius r centered at 
X € and, for x, y G M'^, let Tx^y denote the time when the ball B{y, 1) is 
fully infected in a one-type process with rate A started with Sq = B{x, I). 
Also, write n = (n,0,... ,0). If the radius distribution F satishes ([T]), we 
have 
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lim 

n^oo 


7b,r 


n 


= X a.s. and in L^, 


(3) 


and the same limit is obtained for the time when the single point n is 
infected; see Deijfen et al (2003) or note that, with the modified construction 
of the process, ([3]) follows immediately from Liggett’s subadditive ergodic 
theorem. The reason for considering the time when the entire unit ball 
around the target point is infected (rather than the time at which the point 
itself is infected) is that, in combination with the particular way the model 
is constructed from the Poisson process, it gives rise to the following very 
useful subadditivity property. 


Lemma 2.1 For any x,y,z G 


we have 


T < T 


+ T, 


2,y 


(4) 


Proof: For t > 0, let St denote the infected region at time t starting with 
So = and let St denote the infected region at time t starting with 

Sq = B{z^ 1). We need to show that 


5(y, 1) C (5) 

By definition of we have that B{y, 1) C St^ ^. Hence it suffices to show 
that, for any f > 0 , we have 

St C . (6) 

For t = 0, ([ 6 |) is immediate from the definition of Tx^z- It remains to show 
that the set inclusion is preserved as t increases. The only t’s for which it 
could possibly stop holding are at the times of outbursts in the 5-process. 
Consider the first such outburst at which ([ 6 ]) ceases to hold, denote the 
time (in the 5-process) at which this outburst happens by t', and denote its 
location by x. Since ([ 6 ]) holds for all t < t', we have that at time Tx^z + t' 
the point x has been infected in the 5-process for at least as long as it 
has been infected in the 5-process at time t'. It therefore follows from the 
construction of the process that the 5-process by time Tx^z + 1' has already 
had an outbreak at the point x (or has such an outbreak precisely at time 
Tx,z+t') of the same radius as the outbreak at x in the 5-process. Hence, ([ 6 ]) 
never ceases to hold and consequently holds for all t. In particular, taking 
t = Tz^y yields ([5]), and we are done. □ 

Using the subadditivity in Lemma 12.11 let us prove the following lemma 
which will be crucial in the proof of Theorem 11.11 
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Lemma 2.2 We have 


limsupEfTn _m - To-m] > nX V- (7) 

m^oo 

Proof: Note first that, since E[rn _m] = E[To,n+m] and E[To _m] = E[ro,m], 
we have 

E[Tn__ni — To,_m] = E[ro,n+m “ 7o,m] • 

Hence we are done if we can show that 

limsupE[To,n+m - Eo,m] > nA“V- (8) 

m—>-oo 

To do this, note that, for any positive integer k, E[To,fcn] can be rewritten 
as the telescoping sum 

k 

E[7o,A:n] = ^E[2o,m “ 7o,(i-l)n]- (9) 

i=l 

Suppose now for contradiction that ([8]) fails. Then there exists an e > 0 
such that E[To,in — — e for all but at most finitely many 

positive integers i. Using Q, this implies that 

limsup ^ < nX~^fi — s, (10) 

k—>-oo ^ 

but according to ([3|) we have 

lini = „A- V- (11) 

k^oo k 

Comparing (IlOp and (|lip yields the desired contradiction, and the lemma is 
proved. □ 

3 Proof of Theorem 11.11 

Eirst we quote a proposition that states that the initial configuration is 
basically irrelevant for the possibility of mutual unbounded growth in the 
two-type continuum model. Here P^lrl denotes the probability law of a 
process started with Sq = Ti and 5 q = r 2 . 
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Proposition 3.1 (Deijfen et al 2003 ) Let (ri,r 2 ) and he two 

pairs of disjoint, hounded subsets of with strictly positive Lehesgue mea¬ 
sures. Furthermore, suppose that the radius distribution F has unbounded 
support and satisfies -F(e) >0 for all e > 0. Then 

> 0 ^ (^) > 0- 

Remark 3.1 Proposition 13.11 extends to the case with bounded support as 
well, provided the following (obviously necessary) condition on (P^ j P^) • If 
the radius distribution is bounded by r, then neither of P(^ or P 2 may contain 
an impenetrable layer of thickness r around the other. □ 

Now, to prove Theorem 11.11 we will show that infinite coexistence is indeed 
possible for two infections with the same intensity if the initial sets are 
located sufficiently far away from each other; see Theorem 13.11 In view 
of the above proposition, clearly this implies Theorem 11.11 To simplify 
notation, write 

Theorem 3.1 If F satisfies and Ai = A 2 , then > 0 for large 

n. 

Proof: By time-scaling, we may without loss of generality assume that Ai = 
A 2 = 1. Hence, in what follows, only unit rate processes are considered. To 
avoid superfluous use of superscripts, write = Po,n- 

First, fix e > 0 and pick n large so that 

(i) E[To,n] < (1 + £)npL] 

(ii) Po,n(Po,n < (1 - £)niT) < e, 

which is possible because of ([3|). The symbol ’E’ will throughout the proof 
be used to denote expected value with respect to the marked Poisson pro¬ 
cess underlying the coupling construction described in Section 2. Now, by 
Lemma l2.21 there exist arbitrarily large m such that 

E[rn,-m - Po,-m] > (1 “ (12) 

We will show that, if Po,n(^) = 0 and e is small, then (jl2p fails for large m, 
implying that we must have Po,n(^) > 0. 

Denote by a ball with radius k centered at the point n/2 and, for a 
two-type process started with Sq = B{0, 1) and = B{ii, 1), let 

A\, = {the type i infection reaches points outside B^}. 
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Also, define = A\ H that is, A^ is the event that both infection types 
reach points outside B^. Note that, since A = and A^. D we 

have 

-Po,n(A) = lim Po,n(Afc). 
k—yoo 

Now assume that id),n(A) = 0 and pick m large enough so that -Po,n(Am) < 
1/2. Trivially, 


E[Tn,-m — 7o,--m] — E ( Tn,_m “ 2o,_m )] 

+ E ~ Eo^_m)] ; 

where Ij.} denotes the indicator function. On -'A^, we have T^-m — 
To-m < 0, and hence 


EfTn — Eo,_m] < E (Tn,_m “ Eq _m)] 

— ^ 

where the last inequality follows from Lemma l2.II To bound the right hand 
side from above, write 

E [l{UL}^o,n] = E[To,n] - E • (14) 

By the choice of n, we have E[ro,n] < (1 + e)u^, and so it remains to find 
a lower bound for E[lj_,^^}To,n]- To this end, remember that, also by the 
choice of n, we have Eb,n(Eo,n < (1 — s)nfi) < e. Hence 

E [l{^AL}^o,n] > Eo,n O {To,n > (1 “ e)nii}) ■ (1 - e)nii 

> [Po,nhAl^) - e] • (1 - e)n/r. 

By symmetry, Po,n(-'A)„) = Po,n(-'A^) = Po,n(-'Am)/2, and, using the fact 
that Po,n(Am) < 1/2, it follows that Po,n(“'A)„) > 1/4. Thus 

E [ipALI^o.n] > Q - (1 - e)nfi, 

which in conjunction with p3p and (|14l) yields the estimate 
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E[rn^_in — To^-m] < E[To,n] “ E [l|_,^^|To,n] 

< (1+e)n//- Q-e^ (1-e)n^ 

- (i + 

Now, if e is small, clearly (3/4 +3e) < (1 —2e). Hence, to sum up, assuming 
that Po,n(^) = 0, we have showed that, if e is small and m is large, then 
E[rn,_m—Eo,_m] < (1—2e)n/i. But this contradicts the fact that (fT^ should 
hold for arbitrarily large m and consequently we must have Po,n(^) > 0. □ 
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